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We report results on the static quark potential in two-flavor full QCD. The calculation is performed for three 
values of lattice spacing a -1 « 0.9, 1.3 and 2.5 GeV on 12 3 x24, 16 3 x32 and 24 3 x 48 lattices respectively, at sea 
quark masses corresponding to m^jm p ~ 0.8-0.6. An RG-improved gauge action and a tadpole-improved SW 
clover quark action are employed. We discuss scaling of m p /y/rj and effects of dynamical quarks on the potential. 



1. Introduction 

Interest in the static quark potential calculated 
in full QCD is multi-fold, the foremost being the 
expectation that a flattening of the potential due 
to quark pair creation would be observed at large 
separations. An increase in the Coulomb coef- 
ficient of the potential at small separations is 
another effect expected from sea quarks. A num- 
ber of studies have recently been reported explor- 
ing these aspects of the full QCD static potential 

si- 

Recently the CP-PACS Collaboration has 
started a systematic effort toward a full QCD sim- 
ulation with two flavors of dynamical quarks Q. 
In this article we present results on the static po- 
tential calculated as a part of this program. 

2. Simulations and measurements 

The CP-PACS simulation of full QCD is car- 
ried out for three values of lattice spacing in the 
range a" 1 w 0.9-2.5 GeV and four different sea 
quark masses corresponding to m^/rrip « 0.8- 
0.6. An RG-improved gauge action and a tadpole- 
improved SW clover action are employed to re- 
duce discretization errors. Light hadron masses 
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Table 1 

Parameters of potential measurements. Statistics 
for msL is given in the square brackets. Lat- 
tice scale is from m p — 0.768 GeV at the physical 
point. 



latti 



conf 



12 ;i x24 
0=1.80 
csw = L60 



a. 1 =0.917(10) GeV 



0.1409 
0.1430 
0.1445 
0.1464 



985 
845 
1063 
680 



16^x32 
0=1.95 
csw = 1.53 

1.288(15) GeV 



0.806(1 
0.753(1 
0.696(2 
0.548(4 



0.1375 
0.1390 
0.1400 
0.1410 



"24^x48 
/3 = 2.20 
csw = lA4 



659 [375] 
690 [199] 
659 [169] 
492 [169] 



0.805(1 
0.751(1 
0.688(1 
0.586(3 



= 2.45(9) GeV 



0.1351 
0.1358 
0.1363 
0.1368 



250 
269 
322 
253 



0.800(2 
0.752(3 
0.702(3 
0.637(6 



are calculated at every fifth trajectories g, while 
the static potential calculation uses a subset of 
the configurations as listed in Table |l| where other 
relevant run parameters are also given. 

We extract the static potential from smeared 
Wilson loops by an exponential fit over the range 
T « 0.4-0.8 fm, above which noise quickly dom- 
inates over the signal. The string tension a 
and the Coulomb coefficient a are determined 
by fitting the potential result to a form V(R) = 



• m p (K„ >l =K s J/0.768GeV 

♦ Va/0.440GeV 
i 0.5fm/r 
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A 12 x24, (5=1.8 
■ 16 3 x32, B=1.95 
♦ 24 3 x48, 8=2.2 
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Figure 1. The lattice spacings at /3 = 2.2 deter- 
mined from nip, a and ro as a function of m\a} 
together with quadratic chiral extrapolations. 

Vq — a/R + a ■ R. Errors are estimated by the 
jack-knife method with a bin size corresponding 
to 50 HMC trajectories. 

We also carry out a measurement of the poten- 
tial with 100 configurations in quenched QCD at 
(3 = 2.1 to compare with full QCD results obtained 
at 13 = 2.2. 

At /3 = 1.95, we calculate the static-light meson 
propagator using an exponential smearing func- 
tion for the source. The mass of the static-light 
meson is extracted by making fits over an interval 
T ps 0.3-1.1 fin. 



3. Scaling of m p /y/a 

In Fig. |l] the lattice spacings determined from 
vrip, the string tension a (where we employ the 
value y/a = 0.44 GeV) and the Sommer scale ro 
at each sea quark mass on the 24 3 x 48 lattice are 
plotted as a function of m 2 a 2 . The values extrap- 
olated quadratically in m^a 2 to the chiral limit 
are also shown. We observe that the three lat- 
tice spacings converge within 10% in the chiral 
limit. The consistency of scale observed here at 
a ps 0.08 fm holds also at larger lattice spacings. 

In Fig. H we plot the ratio m p j yfa in the chi- 



ral limit as a function of a p fixed by m p . With 
our choice of the clover coefficient for the quark 
action, the leading 0(g 2 a) scaling violation is ex- 
pected to be small |y. The flat behavior of the 



a p [fm] 

fa in the chiral limit as a 
An estimate 

of the continuum value from a constant fit is also 
shown. 



Figure 2. Ratio m f 
function of a p determined from m p . 



ratio supports this expectation, and also suggests 
that the remaining higher order contributions are 
small as well. We then adopt a constant fit, and 
find m P l \fu = 1.66(4) as an estimate of the con- 
tinuum value. If we employ m p = 0.768 GeV, this 
leads to \[o = 0.463(11) GeV, which is now a pre- 
diction of our full QCD simulation. This value is 
about 5% larger than 0.44 GeV used in Fig. [l| 

4. String breaking effect 

In Fig. |]the potential data V(R) are compared 
with twice the static-light meson mass 2ms l on 
the 16 3 x32 lattice at the lightest sea quark mass. 
The potential result crosses 2msL at around 1 fm 
and continues to increase with R. We find a simi- 
lar behavior at other sea quark masses and lattice 
spacings. Thus our results do not show any in- 
dication of a flattening of the potential expected 
from a pair creation of quarks, similar to those of 
previous studies ||. 

A possible interpretation of this result is that 
the overlap of the Wilson loop operator with 
the state of a static-light meson pair is small. 
In a simplified picture that the contribution of 
the meson pair to the Wilson loop is given by a 
static-light meson propagating along the loop, the 
overlap is given by a factor exp(-EsLR) where 
Esl — msL — fimsL with Stusl ~ Vq/2 is the 
binding energy of meson. Hence the contribu- 
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Figure 3. The potential data on 16 3 x32 lat- 
tice at m w /m p ~ 0.6. The solid and dotted lines 
represent the center value of 2ms l and its error. 

tion of the meson pair dominates over that of 
the string state only for sufficiently large T sat- 
isfying exp(-2£ S L(i? + T)) > cxp(-aRT). For 
numerical parameters appropriate for Fig. ||, we 
find this condition to be met only for T 3> 3 fm 
at R w 1.5fm, which is far too large for practi- 
cal potential measurements. Devising an opera- 
tor which does not suffer from the small overlap 
problem is needed to observe the string breaking 
effect. 

5. Effect on the Coulomb coefficient 

In Fig. [| we plot results for the Coulomb co- 
efficient a obtained on the 24 3 x48 lattice as a 
function of sea quark mass using m^/uip (filled 
circles). The right-most open circle is the re- 
sult in quenched QCD at (3 = 2.7 for which 
a- 1 = 2.217(8) GeV fixed by = 0.44 GeV 
matches with the value a" 1 = 2.22(4) GeV of the 
full QCD run. As expected, the full QCD results 
are larger than the quenched value, and they in- 
crease as the sea quark becomes lighter. 

A simple estimate of the magnitude of shift of a 
from quenched to full QCD may be made in the 
following way||: one starts with the quenched 
value, runs it from fj, — a^ 1 down to /i « 500 
MeV using the 2-loop (3 function with Nf = 0, 
and then runs back to a -1 with Nf = 2. The 
result shown by the horizontal line in Fig. is 
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Figure 4. Coulomb coefficient in full QCD(fillcd 
circle) as a function of m T /m p , as compared to 
the quenched value (open circle) at the same lat- 
tice spacing. Solid and dashed lines represent an 
estimate expected for the full QCD value and er- 
ror. 

consistent with the measured value. Dashed lines 
are uncertainties due to the error of the quenched 
value and the choice of the matching scale, which 
is varied over fi = 400-600 MeV. 
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